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ABSTRACT 

GRB 110721 A was observed by the Fermi Gamma-ray Space Telescope using its two instruments, the Large Area 
Telescope (LAT) and the Gamma-ray Burst Monitor (GBM). The burst consisted of one major emission episode 
which lasted for ~24.5 s (in the GBM) and had a peak flux of (5.7 ± 0.2) x 10 5 erg s 1 cm~ 2 . The time-resolved 
emission spectrum is best modeled with a combination of a Band function and a blackbody spectrum. The peak 
energy of the Band component was initially 15 ± 2 MeV, which is the highest value ever detected in a GRB. 
This measurement was made possible by combining GBM/BGO data with LAT Low Energy events to achieve 
continuous 10-100 MeV coverage. The peak energy later decreased as a power law in time with an index of 
— 1 .89 ± 0. 10. The temperature of the blackbody component also decreased, starting from ~80 keV, and the decay 
showed a significant break after ~2 s. The spectrum provides strong constraints on the standard synchrotron model, 
indicating that alternative mechanisms may give rise to the emission at these energies. 

Key words: gamma-ray burst: general - gamma-ray burst: individual (GRB1 10721A) - radiation mechanisms: 
thermal 

Online-only material: color figures 


1. INTRODUCTION 

Although the emission mechanisms active in the prompt 
phase of gamma-ray bursts (GRBs) are still under debate, 
there is much evidence that the photosphere of the relativistic 
outflow plays an important role in the formation of the observed 
spectrum (e.g., Lazzati & Begelman 2010; Ryde et al. 2010; 
Guiriec et al. 2011; Vurm et al. 2011; Giannios 2012; Pe’er 
et al. 2012). Indeed, a strong contribution from the photosphere 
was predicted on physical grounds in early works by Goodman 
(1986) and Paczynski (1986) but this was not considered a viable 
model since the observed spectra are, in general, non-thermal. 
It was, however, realized that the photospheric emission should 
be accompanied by non-thermal emission from optically thin 
regions (e.g., Meszaros et al. 2002), and that the photospheric 
emission can be enhanced and modified from a Planck function 
by energy dissipation at moderate optical depths (Rees & 
Meszaros 2005). 

The Fermi Gamma-ray Space Telescope has made obser- 
vations that support this view. The very bright GRB090902B 
(Abdo et al. 2009) was observed to have a narrow and steep 
spectral component that may best be attributed to the photo- 
sphere (Ryde et al. 2010; Zhang et al. 2011). At later times 
during the burst this spectral component widened into a broader 
feature. This implies that the spectrum of photospheric emission 
must be able to have a variety of shapes, not only a Planck func- 
tion (Ryde et al. 2011). A probable explanation for this is the 
existence of sub-photospheric energy dissipation. Indeed, some 
models predict that a thermal peak could be well described by 
the Band component (e.g., Beloborodov 2010). Furthermore, 
many bursts have shown signs of a subdominant photospheric 
component. For instance, GRB090820A exhibited two spectral 
peaks: a peak related to a blackbody spectrum with a tempera- 
ture of ~40 keV, and a peak modeled by a Band function (Band 
et al. 1993) at E p ~ 1 MeV (Burgess et al. 2011). 

Only a few spectra having peak energies of a few MeV 
have previously been reported. Gonzalez et al. (2009) presented 
GRBs observed over the energy range 0.02-200 MeV, through 
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a joint analysis of data from BATSE and the EGRET TASC on 
the Compton Gamma Ray Observatory ( CGRO ). The maximal 
value they found was 5 ± 2 MeV in GRB981203. Similarly, the 
PHEBUS experiment on board Granat, with an energy range of 
0. 1-100 MeV, identified spectral breaks in other bursts at around 
2.4 MeV in addition to the low-energy break (Barat et al. 2000). 
Fermi observations of GRB090510 showed a peak energy of 
~5 MeV (Ackermann et al. 2010). 

In this Letter, we study GRB110721A, which besides being 
very bright during the first ~8 s (a few x 10~ 5 erg s -1 cm -2 ) 
has an initial peak energy of a record breaking 15 MeV — 
observations made possible by Fermi’s exceptional spectro- 
scopic capability. Moreover, the y- ray emission is dominated 
by a FRED (fast rise, exponential decay) pulse, which makes 
this burst ideal to compare with previous studies of such pulses. 
Apart from the Band component, we identify, with high signifi- 
cance, an additional component which we model using a Planck 
function. We interpret this as a photospheric component. 

2. OBSERVATIONS 

On 2011 July 21, the Gamma-Ray Burst Monitor (GBM; 
Meegan et al. 2009) and the Large Area Telescope (LAT; 
Atwood et al. 2009) on board Fermi detected high-energy emis- 
sion from GRB1 10721A (GCNs 12187, 12188). The burst po- 
sition was triangulated by the Interplanetary Network (IPN), 
which returned an error box centered in (R.A., deck) = 
332?46, -38:63 (J2000) and approximately 1°3 x 0:4 wide 
(GCN 12195). The intersection of the IPN box with the LAT 
error circle gives an area of ~ 1 200 arcmin 2 centered on 
(R.A., deck) = 333°2, — 38°5. We adopted the latter posi- 
tion in our analysis. Figure 1 shows a composite count light 
curve from the various detectors on Fermi : the Nal detectors 
(8-900 keV), the BGO detector (200 keV-40 MeV), and the 
LAT (P7V6_Transient class events; Atwood et al. 2009). The 
most energetic photon was detected at 4.50 s, had an energy 
of E = 6.3 ± 0.6 GeV, and was associated with the GRB 
with a high (>0.9 using the gtsrcprob Fermi Science tool 54 ) 
probability. The light curve of the burst is dominated by a sin- 
gle and exceptionally bright FRED-like emission episode with 
790 = 24.5 s, the time during which 90% of the emission is 

54 http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
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Figure 1 . Composite light curve of GRB1 10721 A. (a)-(d) Count light curves from different energy ranges (Nal 6, 7, 9, 11, and BGO 1 detectors), (e) Individual LAT 
photons. Time intervals are indicated by green lines; the red line shows the trigger time. Filled circles in (e) indicate >90% probability of association with the GRB. 
The time is relative to the GBM trigger. 

(A color version of this figure is available in the online journal.) 


received. The peak of the energy flux over the energy range 
8 keV-1 GeV occurs at 0.3 s relative to the GBM trigger with 
(5.7 ± 0.2) x 10 -5 erg s _1 cm -2 . 

In Figure 1 we also present the light curve of the LAT 
Low Energy events (LLE, 30-130 MeV). The LLE data are 
produced from a non-standard LAT analysis which for bright 
sources provides large effective area at low energies, joining the 
LAT and GBM energy ranges (Pelassa 2011). An observation- 
specific response matrix is generated for each time range of data 
using a Monte Carlo simulation of the LAT, which increases 
the spectral capabilities of the Fermi data. We note that the 
LLE light curve for this burst is peculiar since it peaks before, 
and its duration is significantly shorter than, the GBM light 
curve, in contrast to what is typically observed in other bursts 
(Fermi LAT Collaboration, in preparation). Finally, in Figure 1 
we explicitly show the Nal light curve in the narrow energy 
range of 8-100 keV. Interestingly, at these energies the light 
curve differs in that a second peak appears at ~2 s. These facts 
will be discussed further in Section 4. 

We have performed a standard analysis with the RMfit 4.0 
package (Mallozzi et al. 2005) using Time Tagged Event (TTE) 
data from the Nal 6, 7, 9, 11, and the BGO 1 detectors. We also 
include the LLE and LAT data. 

3. SPECTRAL BEHAVIOR 

A Band function fits the integrated emission poorly during 
the pulse (—0.32 to 8.38 s). The Castor C-statistic (C-stat; 
Arnaud et al. 2011) for the fit is 1078 for 618 degrees of freedom 


(dof). Based on the earlier detection of a blackbody component 
in other Fermi bursts (Ryde et al. 2010; Guiriec et al. 2011; 
Burgess et al. 2011), we test whether the fit can be improved 
by including a blackbody. The fit improves significantly to 
C-stat/dof = 901/616 (corresponding to >5 ct significance). 
Compared to the Band-only fit, the peak energy of the time- 
integrated Band function shifts from E p = 1120 ± 60 keV to 
2400 ± 170 keV. Since the blackbody component is expected 
to evolve (see analysis below and Ryde & Pe’er’s 2009) its 
time-integrated spectrum is expected to be better characterized 
by a multicolor blackbody (mBB; Pe’er & Ryde 2011). Using 
such a component instead of the blackbody yields a significant 
improvement of the fit C-stat/dof = 871/615. Our preferred fit 
for the time-integrated spectrum is thus the Band + mBB model 
as shown in Figure 2. 

We also analyze the time-resolved spectra using data up to 
130 MeV (i.e., the LLE data). In order to study the spectral evo- 
lution we need a temporal resolution as high as possible. How- 
ever, since the signal-to-noise ratio (S/N) then also decreases, 
there is a trade off in our choice of binning. We require an S/N 
of 40 in the most strongly illuminated GBM detector (Nal 9) 
in the energy range 8.0-100 keV. This ensures that the spectral 
fits are well constrained over the energy range where the black- 
body component was found in the time-integrated spectrum. 
This procedure divides the bursts into eight time bins which are 
indicated in Table 1 . 

As for the time-integrated spectrum, we find that the addition 
of a blackbody to the Band component improves the fit. In 
the first seven bins the C-stat value decreases by a significant 
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Table 1 

Fit Results Using a Band Function Only and a Band+Blackbody Model 




Band Function Only 




B and+B lackbody 



Time (s) 

Fpeak (keV) 

a 

P 

C-stat/dof 

kT (keV) 

Speak (keV) 

a 

P 

C-stat/dof 

-0.32-0.83 

0.83-1.54 
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Figure 2. Spectral fit and residuals of the time integrated emission (—0.32 to 
8.38 s) and the best-fit model, Band + mBB. The top panel shows a vF v spectrum, 
the middle panel a count spectrum, and the bottom panel the residuals of the fit. 
(A color version of this figure is available in the online journal.) 


amount (see Table 1 ). Using an assumed Band function spectrum 
as the null hypothesis, and Band+Blackbody as the alternative 
hypothesis, we determine the significance of the additional 
blackbody component. The distribution of the test statistic AS 
was investigated through Monte Carlo simulations. We find that 
an improvement of AC-stat = 30 when adding a blackbody 
component to the fit, i.e., adding 2 dof, corresponds to a 1 0‘ 7 
probability of not having a real blackbody component in the 
spectrum. Similarly, AC-stat = 20 corresponding to a 10™ 
probability and AC-stat = 10 to a 1 0 3 (3a) probability. In 
particular, the normalization of the blackbody component is 
constrained for the first seven bins. We tested the robustness 
of this result by using another time bin selection based on 
Bayesian blocks (Scargle 1998). The results are fully consistent. 


We therefore conclude that an extra component is significantly 
detected in the time-resolved spectra up to ~6 s. 

In Figure 3 we plot the observed temperature of the blackbody 
as a function of time. First, it is apparent that there is a strong 
evolution of the temperature. Second, there is a notable break in 
the decay of the temperature. The energy flux in the blackbody 
component is approximately 5% of the total flux and it peaks at 
around 2 s, when it reaches 10%. In the energy band around the 
peak it is much stronger, and well above the systematic errors 
(10%; Bissaldi et al. 2009). We note that this peak is notably 
different from the peak of the energy flux pulse, which occurs 
at 0.3 s, but coincides with the temperature break. Moreover, 
it coincides with the second peak in the Nal count light curve 
(Figure 1) which is mainly due to photons below 100 keV. This 
suggests that this second peak is associated with the blackbody 
component. The temperature is well fitted with a broken power 
law (see Ryde 2004) in time with a break at 2.3 ± 0.2 s. Note 
that such a break does not appear in the decay of E p (Figure 4). 
The power-law indices before (after) the break are —0.30 ±0.13 
(—1.66 ±0.15). This behavior is quantitatively the same as what 
was found for the photospheric component in CGRO BATSE 
bursts (e.g., Ryde & Pe’er 2009). 

In order to study the evolution of the Band component, we 
allow for a higher temporal resolution since it is the dominating 
spectral component. We adjust the time bins to provide a 
minimum S/N of 30 in the counts spectrum from the Nal 
9 detector using the full energy range (8.0-860 keV), which 
increases the temporal resolution while still maintaining good 
spectral constraints. Since the burst is initially very strong in 
the LLE data and spectral evolution is initially very rapid, we 
divide the first time bin into two. This gives us 13 time bins. 

A striking feature in this burst is the unusually high value 
of the peak energy. The highest value is measured in the first 
time bin (—0.32 to 0.0 s). Here the peak energy is 15 ± 2 MeV, 
the low-energy slope has a photon index a — —0.81 ± 0.08, 
the high-energy power law has /3 = — 3.52° 0 g. The fit has 
C-stat/dof of 679/608. Note that a blackbody component cannot 
be constrained by the data in this time bin, but we include it in 
the following bins. The peak energy of the Band component 
decreases monotonically with time; E p = A p \(t — to) d , where 
d = —1.22 ± 0.13 and t 0 = -0.46 ±0.10 s (Figure 4). 
Attempting to fit the data before and after 2.3 s (the time 
where the temperature evolution shows a break) separately gives 
indices that are compatible within errors. We note that to is 
consistent with the onset of emission in the LLE range. 

After a few seconds E p reaches a few hundred keV, typical 
for other GRBs (Kaneko et al. 2006). Consistent results are 
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Figure 3. Blackbody temperature, kT, decays as a broken power law (fit 
function from Ryde 2004). Circles correspond to the binning based on the data 
below 100 keV, filled/open indicating 5a /3a significance of the blackbody 
component. Gray squares show results from the higher time resolution binning. 
(A color version of this figure is available in the online journal.) 

obtained when using other temporal binnings. Figures 3 and 4 
show the values found for both the coarser (red points) and finer 
(gray squares) time resolution. The two sets of points are clearly 
consistent. 

The binning used for the time-resolved analysis results in very 
few data above 100 MeV in each bin. We therefore also studied 
the time-integrated spectra in three broad time intervals in order 
to increase the signal at the highest energies: [—0.325. 1.5] s, 
[1.5,5] s, [5,20] s. For each of these bins an extrapolation 
of the best-fit (Band+BB) model of the data below 100 MeV 
is consistent with the LAT data. The high-energy emission 
continues for more than 200 s after the burst trigger (Tgo for 
the LAT emission is 253 s) and decays as a power law with 
index —0.95 ± 0.04, a typical value for the temporally extended 
emission observed in other bursts ( Fermi LAT Collaboration, in 
preparation). 

4. DISCUSSION 

The value of E p = 15 ± 1.7 MeV found in the initial time bin 

in GRB 110721 A is the highest measured in a GRB spectrum. 

Capturing such a high initial peak energy is a testament to 

the importance of acquiring good quality BGO and LLE data 
in bursts; this Letter provides the first such realization of this 
contribution to the GRB paradigm. 

Bursts with hard spectra below their vF v peaks are frequently 
observed. Those with indices a > —1.5 below this peak 
cannot possess electrons that radiate synchrotron emission in 
the expected fast-cooling regime, within this spectral window; 
this is the so-called fast-cooling a index limit (Preece et al. 

1998). GRB 1 10721 A is just such an example, and its spectrum 
can be consistent with optically thin synchrotron emission 
only if the non-thermal electrons are cooled on timescales 
much longer than those on which they are injected/replenished. 
Demanding that the cooling break should lie well above 15 MeV, 
combined with the requirement that the observed vF v peak at 

1 5 MeV in the first time window corresponds to the synchrotron 
peak, constrains the emission region to be smaller than typical 
photospheric radius scales (F. Ryde et al., in preparation), i.e., 

10 1() -10 4 * * * * * * * * 13 cm when accounting for GRB bulk motions (e.g., 

see Rees & Meszaros 2005). This is a significant restriction 
on synchrotron emission models that invoke dissipation zones 
outside the photosphere. Therefore, it is desirable to entertain 
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Figure 4. Evolution of the E p of the Band function component. Symbols as in 
Figure 3. 

(A color version of this figure is available in the online journal.) 

alternative explanations for the very “blue” emission seen in 
GRB 1 10721 A, e.g., Compton scattering of thermal photons. 

The cleanest signature of the emission in GRBs is found 
in smooth FRED emission pulses (Ford et al. 1995). Ryde 
(2004) identified a few such pulses in the CGRO BATSE catalog 
which were consistent with having a Planck function spectrum 
throughout their durations over the observed spectral range 
25-2000 keV, suggesting a photospheric origin of the emission. 
It was further argued in Ryde (2005) that the spectral break 
in apparently non-thermal spectra can be interpreted as the 
photospheric emission peak. While the peak is modeled by a 
Planck function, an additional power-law component is needed 
to account for the non-thermal character of the spectrum. More 
importantly, the temperature of the blackbody in all pulses is 
observed to have a characteristic recurring behavior: it decreases 
following a broken power law in time (see further Ryde & 
Pe’er 2009). These results suggest that photospheric evolution 
in GRBs may exhibit well-defined characteristics. 

The evolution of kT seen here is similar to what was 
observed by BATSE over individual GRB pulses where a 
thermal emission component was identified (Ryde 2004; Ryde 
& Pe’er 2009). Moreover, if this burst had been observed by 
BATSE it would be best modeled by a blackbody and a power 
law, which is apparent by limiting the energy range to that of 
BATSE. GRB 110721 A therefore confirms the BATSE results 
on the behavior of the photospheric emission. 

Finally, the two peculiarities of the 30-100 MeV emission, its 
early onset and short duration, can be explained by the fact that 
the emission originates in the same component as detected by the 
GBM. This component initially dominates the LLE range due 
to the exceptionally high value of E p but moves into the GBM 
range as the burst evolves. Therefore, the emission detected 
in the LLE range precedes, and has shorter duration than, the 
emission seen at lower energies; this mirrors the behavior seen in 
BATSE bursts with FRED-like light curves (Norris et al. 1996). 

In summary, for GRB 11 0721 A we have shown that 
(1) initially the peak energy has an extreme value E p = 15 MeV, 
the identification of which was made possible by the high- 
quality LLE data. This value combined with the fact that 
a > —1.5 cannot be explained by the standard optically thin 
synchrotron model. (2) In addition to the Band function a nar- 
row spectral component (consistent with a blackbody) is statis- 
tically required by the data. (3) The spectral evolution of this 
blackbody is similar to what was found in BATSE pulses hav- 
ing a photospheric component. (4) The blackbody component 
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evolves differently compared to the Band component. (5) The 
second peak in the Nal count light curve coincides with the 
time when the blackbody component flux is the highest. These 
facts provide strong evidence for the existence of photospheric 
emission in GRB1 10721 A. 
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